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It has been made an extention of the conventional theory based on the assumption of the well isolated Zhang- 
.Rice singlet to be a first electron-removal state in dielectric copper oxide. One assumes the photohole has been 
"localised on either small (pseudo)Jahn- Teller polaron or large non-adiabatic polaron enclosed one or four to five 
• CuO^ centers, respectively, with active one-center valent C'Aig — ^'^Eu) manifold. In the framework of the cluster 
'model we have performed a model microscopic calculation of the k-dependence of the matrix element effects and 
.photon polarization effects for the angle-resolved photoemission in dielectric cuprate like Sr2Cu02Cl2- We show 
'that effects like the "remnant Fermi surface" detected in ARPES experiment for Ca2Cu02Cl2 may be, in fact, 
,a reflection of the matrix element effects, not a reflection of the original band-structure Fermi surface, or the 
"strong antiferromagnetic correlations. The measured dispersion-like features in the low-energy part of the ARPES 
.spectra may be a manifestation of the complex momentum-dependent spectral line-shape of the large PJT polaron 
"response, not the dispersion of the well- isolated Zhang-Rice singlet in antiferromagnetic matrix. 

Introduction 

Angle-resolved photoemission spectroscopy (ARPES) is considered to be a key experiment to elucidate a number 
■ of principal issues of electronic theory related to the unconventional properties of cuprates, manganites, nickellates, 
'bismuthates, and other strongly correlated oxides. These are the quasiparticle dispersion, the Fermi surface, 
I (pseudo)gap behaviour and other ones. Many researchers consider namely ARPES data to be a main evidence in 
' favour of either model of electronic structure and energy spectrum. 

However, the current practice of ARPES activity in oxides does not justify these expectations. Up to now the 
ARPES measurements are treated as a rule in the framework of oversimplified " three-stage" approach. The first 
."experimental" one relates straightforwardly to measurements of the {E, k)-dependencies of photocurrcnt intensity. 
"Any problems here are related to a sample preparation, experimental setup, and attempts to avoid such external 
.parasitic phenomena as a charging effect. The very "dangerous" second "experimental-theoretical" stage implies the 
' preliminary treatment and often incorporates a " hidden" interpretation based on a number of seemingly inviolable 
.statements like: 

1. ARPES data in a wide spectral range could be self-consistently described by the conventional band models. 
Electronic structure and energy spectrum near the Fermi level are formed by Landau quasiparticles with Fermi 
'statistics which fill an energy band. 

2. The peaks position for the photocurrcnt intensity links to the peaks in the quasiparticle spectral density. In 
other words, ARPES essentially measures the one-particle spectral function of the initial state. 

3. Some authors go further and make use of the standard method which provides the information on the Fermi 
'surface from ARPES data in the conventional band metals. First, one find the occupation probability, 7i(k), by 

integrating the ARPES spectral function A(k, i?) over energy. Experimentally, one choses an energy window for 
'integration, thus becoming the relative rt(k). Then, the drop of the relative n(k) is used to determine the Fermi 
,. surface. In practice, one defines fci? as the locus of points of maximum gradient of rt^r or as the point of steepest 
'descent in the relative n(k). Naturally, this method implies the simple metallic- like electronic structure, where the 

identification of a Fermi surface is convincing. 

In our opinion, this "band-like" paradigm supported only by the simplest model like that of free electrons, turns 
out to be at least questionable in all points in the case of strongly correlated oxides. Firstly, an ability of band 
'models, even modified like LDA-I-U [Q, to yield a relevant description of electron spectra for strongly correlated 
oxides is merely a misleading, at present there are not convincing examples for such a description, in particular, 
for the lowest-lying occupied band, or the so-called first electron-removal states. The typical errors in locating 
the final states for ab initio band structures are generally expected to be of the order of a few eV's (!). Large 
number of experimental data evidences in favour of unconventional "non-Landau " nature of quasiparticles and 
occurrence of unusual correlations in cuprates and other oxides. Traditional interpretation of ARPES data like a 
simple accordance between the photointensity peak position and the quasiparticle spectral density peak implies a 
full neglect of several factors each of which is capable to result in crucial reconsideration of the ARPES data: i) the 
matrix element effects, or the intensity dispersion; ii) surface effects; iii) effects of multiple scattering; iv) effects of 
finite lifetimes of the initial and final states; v) effects of coupling with phonons, spin excitations and other possible 
degrees of freedom; vi) effects of configurational interaction; vii) effects of electron inhomogeneity, for instance, in 
doped cuprates. For insulating samples it is difficult to escape the charging effects. 

Importance of matrix element effects and the ARPES intensity dispersion for different cuprates was underlined 
earlier and a role played by configurational interaction was partially illustrated in Ref. 1^ for Sr2Cu02Cl2- 



The hole polaron formation due to electron-phonon coupHng and its manifestation in ARPES was demonstrated 
by Alexandrov and Dent |^ for cuprates like Y123 and Y124. One should be noted that all the above mentioned 
citations relate to the nontraditional approach in ARPES. 

The third, purely theoretical stage implies a qualitative and quantitative description of the dispersion law and 
other quasiparticle properties obtained on the foregoing stage and to be considered as "experimental data" (?!). 
Unfortunately, too often this stage looks like a formal fitting with either model theory, and with no critical inspection 
of experimental data. Here one might mention, for instance, numerous papers with theoretical treatment of 
"experimental" ARPES data for Sr2Cu02Cl2 in the framework of various extended t — J, t~t' — J, t — t' — t" — J 
models §,|lO|. 

Summarizing, we see this "three-stage" approach could result in a natural doubts as to conclusions based on 
such an oversimplified interpretation. 

As a convenient demonstrative model system for ARPES in the insulating layered copper oxides one might be 
chosen oxychlorides Sr2Cu02Cl2 and Ca2Cu02Cl2, which are isostructural to famous 214 system La2CuOi. The 
oxychloride Sr2Cu02Cl2 is one of the most popular model system for insulating phase of the high-Tc cuprates and 
is intensively studied both experimentally and theoretically. 

In this tetragonal antiferromagnetic with nearly ideal Cu02 planes there are chlorine atoms instead of apex 
oxygens with considerably larger Cu — Clapex separation (2.86A) than that of Cu — Oapex (2.42A) in La2CuOi. 
Hence, in Sr2Cu02Cl2 one appears a real opportunity to examine the Cu02 plane states, both copper and oxygen, 
without "parasitic" contribution of apex oxygens. At present there are a rather large number of experimental 
data for Sr2Cu02Cl2 obtained with the help of optical spectroscopy X-ray photoemission (XPS) |12|, ultra- 
violet photoemission (UPS) X-ray absorption spectroscopy (XAS) 0,^5), electron energy loss spectroscopy 
(EELS) [l^-^l^ , angle-resolved photoemission spectroscopy (|J|,^,^,|2| . Similar experimental ARPES results were 
obtained recently by F. Ronning et al. for Ca2Cu02Cl2 which is a full analog of Sr2Cu02Cl2- 

The authors |^ state: "This is a measurement of the dispersion of a single hole in an antiferromagnetic back- 
ground, a problem that has been heavily investigated theoretically" . It seems, this is a very strong and unambiguous 
statement, which oversimplifies the problem and does not account of many alternative scenarios. Moreover, the 
experimental data |^ could hardly considered to be a reliable basis set for any decisive conclusions. Indeed, the 
spectral dependence of the photocurrent intensity shows a complex nature of the low-energy photo-hole states with 
a strong and dispersive contribution of the high-lying states. 

It is very important to emphasize that, while the ARPES spectra in various oxides displayed as the photointensity 
against i?, k curves show the dispersion-like features (see, e.g. Fig.l, Ref. [^, or Fig. 3, Ref. (20)), the "dispersing" 
states which " peak positions" are plotted are extremely broad, with width comparable to binding energy, and these 
simply cannot be thought of as quasiparticles. This general point is true at all k's. In addition, in the most part of 
the Brillouin zone there are no well-defined ARPES peaks at all. Moreover, the observed peak position may rather 
strongly depend on the energy resolution. 

The similar to Sr2Cu02Cl2 experimental results were obtained by F. Ronning et al. for Ca2Cu02Cl2- 
The improved spectral quality allows authors to reveal a steep drop in spectral intensity across a contour that is 
close to the Fermi surface predicted by the band calculation. They concluded that the Fermi surface, which is 
destroyed by the strong Coulomb interactions, left a remnant in this insulator with a volume and shape similar to 
what one expects if the strong electron correlation in this system is turned off. The lowest energy peak exhibits 
a dispersion with approximately the \coskxa — coskya\ form along this remnant Fermi surface, in other words 
the strong correlation effect deforms this otherwise iso-energetic contour (the non-interacting Fermi surface) into 
the form that resembles the d-wave like pseudogap dispersion with a very high energy scale of 320 meV [20| . 
The authors j20| consider the d-wave like dispersion of the parent insulator to be the underlying reason for the 
pseudogap in the underdoped cuprates. They follow straightforwardly the simplest model approach to ARPES 
and electronic structure of the photo-hole, and do not consider either alternative approaches, albeit they mention 
some "parasitic" effects of matrix elements and photon polarization, especially, if these effects hardly keep within 
their generic model. The approach |Q is demonstrative for many papers on ARPES, in particular, addressing the 
Fermi surface problem. 

At present there are several papers with experimental ARPES data on Sr2Cu02Cl2 and Ca2Cu02Cl2 systems 
, ^ , p2| , p0[ distinguished by different energy resolution from 25 meV to 105-115 meV different photon 
energy from 10 eV to 80 eV [pPp^, light polarization pl] , p3| |, and the measurement temperature. Despite the 
many common features there are some important departures. So, contrary to findings |9j, the strong "quasiparticle" 
dispersion was observed in Sr2Cu02Cl2 for the (tt, 0) direction with strong and unusual dependence of the peak 
amplitude on the light polarization |^l|. Kim et al. |^2| performed the ARPES measurements in Sr-oxychloride 
at 150 K, well below the Neel temperature, and have confirmed that the spectrum in the (tt, tt) direction consists 
of a single relatively sharp peak near (7r/2,7r/2). However, the spectra along (tt, 0) direction are very broad and 
consist of at least two peaks separated by about 0.4 eV . Just recently, using the ARPES spectra of Sr2Cu02Cl2 
as an example, the authors j^,^ have experimentally demonstrated a significant impact of electron-photon matrix 
elements on both the relative spectral intensity and the shape of a low-energy feature in ARPES spectrum. 

Many experimental ARPES data and conclusions on Sr2Cu02Cl2 known up to now were revised very recently 



in a detailed investigation by C. Diirr et al. ||2^. The authors observed the marked oscillating photon-energy 
dependence of the photoemission signal of the first electron-removal state, which was attributed to the diffraction 
of the photoelectron wave on the c-axis periodically arrayed Cu02 planes. They found a strong polarization 
dependence in ARPES spectra along the high-symmetry directions, compatible with that expected for a Zhang- 
Rice singlet. The ratio between the coherent and incoherent spectral weight in ARPES spectra near the first 
electron-removal state appears to be photon-energy dependent. Among several mechanisms of this puzzling effect 
the authors |^3| mention the possible contribution of the low-lying electronic states other than the Zhang-Rice 
singlet. 

One of the remarkable spectral features clearly revealed in the ARPES measurements for both insulating cuprates 
[ppC|] is a significant low-energy spectral weight in the F (0,0) point that accordingly to symmetry of the electric- 
dipole matrix elements implies a significant weight of the purely oxygen odd Cu-symmetry photo-hole state. How- 
ever, in the framework of traditional approach one prefers to take no notice of this feature. 

Our paper has not for an object the elaboration of the general theory of ARPES in the strongly correlated oxides. 
The authors would like to emphasize on the possibly simple but real examples an importance of some factors such 
as matrix element effects, formation of the hole non-adiabatic polaron, often to be disregarded in the framework 
of conventional approaches to interpretation of ARPES data in copper oxides. At the same time the paper could 
be considered as a first step in elaboration of an original theory of ARPES in the strongly correlated oxides. 

The paper is organised as follows. In Section I we present a brief summary of results of quantum chemical 
modelling the electronic structure and energy spectrum for the CuO\~ and CuO\~ centers with one and two holes, 
respectively. One considers the formation of small and large non-adiabatic polarons. It should be mentioned here 
that the concept of polarons led to the discovery of the copper oxide superconductors. Section II incorporates an 
analysis of the generic expression for the photoemission current intensity in the framework of the PJT polaronic 
approach. In Section III we present a model calculation of dipole matrix elements which describe a transition of 
a bound 7/i electron to unbound state. An illustrative modelling of the k-, and polarization dependence for the 
matrix element part of the ARPES intensity is given in Section IV. 

1. Electronic structure of the hole state in cuprates. Prom small to large PJT 
polaron 

One of the peculiar properties of the doped cuprates is a cross-over from localized to itinerant electronic behaviour 
with an inhomogeneous distribution of electronic states which description in a wide range of compositions from a 
single-valent antiferromagnetic insulator to a mixed valent bad metal is a fundamental problem of the solid state 
physics. 

Unfortunately, at present there is no a general consensus concerning the nature of valent electron and hole states 
which form optical and electron (PES, XPS, EELS) spectra even in insulating cuprates hke La2Cu04, YBa2Cu30e, 
Sr2Cu02Cl2- It is becoming increasingly difficult to reconcile experimental results with the expectations of the 
simple band models and Fermi liquid theory. 

One of the firmly established facts relates the Cu02 plane character of the valent states with the hig{dx2-y2) 
hole ground state for the CuO\~ center. Various spectroscopic methods link the low-energy excitations in the 
range E > 2 eV to the charge transfer 02p Cu3d within the Cu02 plane. At a first glance, optical and 
electron spectra manifest both localised and delocalised, or band character of electronic states. Lack of reliable 
interpretation of optical and electron spectra in strongly correlated oxides appears to be a result of difficulties in 
theoretical description of strong covalency and strong intracenter and intercenter correlations. 

The parent cuprates such as CuO, Sr2Cu02Cl2, La2Cu04, YBa2Cu30Q provide typical examples of systems 
with strong local correlations when dielectric antiferromagnetic ground state is mainly specified by potential energy 
of electron-electron coupling. Kinetic energy would prefer formation of half-filled band with typically metallic 
behavior. Occurrence of the dielectric antiferromagnetic phase of the parent cuprates itself and wide opportunities 
of the interacting Cu04-centers model in explanation of many optical and electron spectra for various copper 
oxides could be considered to be a convincing evidence in favor of quasiparticle states localised predominantly 
on the CM04-clusters which do not obey conventional band model description. This is a common place as for 
Cu3d — 02p-hybrid states, but purely oxygen non-bonding states are currently described as the band ones. Such a 
discrimination is partially based on a rather simplified relation between correlation and Slater electrostatic integrals 
like Fo{3d3d) > Fo{2p2p), and typical, for instance, for the Anderson impurity model. 

However, correlation effects appear to be of particular importance for oxygen hole states which was pointed out 
by Hirsch et al. in their theory of "anionic metal". A simple model for localisation of the oxygen holes and 
some other puzzling consequences of conjectured multiplet structure for anionic background were considered 
recently . A strongly correlated behaviour of the 02p-holes to some extent gives them equal rights with 3c?-holes 
and leads to reconsideration of many conventional approaches and models which ignore strong 02p-correlations, 
for instance, the Anderson impurity model which considers purely oxygen states to be band-like. 

So, a quantum-chemical cluster approach appears to be more relevant for description of strongly correlated 
oxides as it allows to account for correlation effects both for cation and anion states by the most optimal manner. 
It seems, the parent copper oxides can be treated by a conventional quantum-chemical ligand-field theory, so a 



localized description can be considered to be a relevant starting point for modelling the electronic structure both 
in ground and excited states. 

First of all, the modelling implies a choice of the effective basis set of atomic orbitals, such as Cuid- and 02p- 
orbitals in the cuprates. Then we introduce the CUO4 center being the minimal atomic cluster that posesses 
the relevant point symmetry D^h and incorporates the main CuSd — 02p covalent bond. This center could be 
considered as an elementary building block to effectively model the electronic structure both in ground and excited 
states. Along with the symmetry and covalency the CUO4 center allows to account for the electronic and vibronic 
correlations. 

1.1. Electronic structure and energy spectrum for the CuO\~ and CuO\~ centers in cuprates 

Slightly distorted tetragonal CuOi center is the only common cell of crystalline and electronic structure in a 
wide variety of high-T^ copper oxides and related parent systems. Restricted atomic basis of five Cu3d and twelve 
02p atomic orbitals for the CuO^ cluster with D^h symmetry yields seventeen symmetrised orbitals with aig, 02^, 
big, b2g, eg (even - gerade) and a2u, f'2u, e„ (odd - ungerade) symmetry. The even Cuid orbitals with aig^id^^), 
big{^dx2-y2), b2g{idxy), eg{3dxz, -idyz) symmetry hybridize with even 02p orbitals with the same symmetry thus 
forming the appropriate bonding 7^ and antibonding 7" molecular orbitals. Among the odd orbitals only e„(fT) 
and (tt) ones hybridize each other due to a strong 02p — 02p coupling to form bonding and antibonding purely 
oxygen molecular orbitals and e", respectively. Purely oxygen a2g,a2u,b2u orbitals are nonbonding. All the 
plane molecular orbitals could be subdivided to a (ai^, e„(cr)) and tt (a2g, 623, eii(7r)) ones, depending on the 
orientation of 02p orbitals. 

Quantum-chemical cluster approach allows the optimal account for electrostatic correlations and description of 
a rather complex electronic structure of the two-hole CuO\~ center. Firstly, one should note an importance of 
configurational interaction, which could be rather simply illustrated by the model cluster calculation ||29|] . Indeed, 
the wave function and energy of the ground state term ^Aig (Zhang-Rice singlet) exhibits an essential coupling for 
three configurations like bigbig-. 

\^Aig Q.mblgf > +0.55|6?,6j, > -QAQWigf > . 

If to consider the ^Aig singlet originating from the only {b\g)'^ configuration and to represent its energy as follows 

ECAig) = 2sib1g) + U + AU, 

with U being the contribution of the hole-hole b'lg — b\g interaction, and AJ7 that of configurational interaction, 
then 

C/ -f AC/ = 4.7 - 3.5 = l.2{eV), 

which value is puzzlingly small as compared with the bare value Ud =^ A + AB + 3C ~ 10 eV for purely atomic 
dl2_y2 configuration. 

A concept of the well isolated Zhang- Rice singlet to be a ground state of the two-hole CuO\^ center, is a guideline 
of many popular model approaches js^]. Namely with ZR singlet one associates the first electron- removal state 
and the lowest-lying features in ARPES spectra. At the same time, a number of experimental data and theoretical 
models evidence a more complicated structure of the valent multiplet for the two-hole CuO\~ center. 

A model of the valent ^ Aig — ^ multiplet, developed in Refs. |p9| , ^ , ^ , implies a quasi-degeneracy in the ground 
state of the two-hole CuO\^ center with two close in energy ^ Aig and^Eu terms of b\g and bigCu configurations, 
respectively. In other words, one implies two near equivalent allocations for the additional hole, either to the 
Cu3d02p hybrid big state, or to purely oxygen e„ state with peculiar Cu^+ — Cu^'^ valence resonance (see Fig.l). 
Occurrence of the localized purely oxygen like states is provided, in particular, by the specific properties of the 
"non-rigid" anionic 02p^ background p^ , p8t . 

The model approach under consideration is based both on microscopic quantum chemical study of the model 
copper-oxygen clusters p9|-^ and a large variety of experimental data. To the best of our knowledge one of the 
first quantitative conclusions on the competitive role of the hybrid copper-oxygen big{dr^2_y-2) orbital and purely 
oxygen 02pt^ orbitals in formation of valent states near the Fermi level in the Cu02 planes has been made by Jiro 
Tanaka et al. (see also more later publication 

In a sense, the valent {b\gYAig — {bigCuYEu manifold for the hole CuO\~ center implies an unconventional 
state with Cu valence resonating between Cu^^ and Cv?^ , or " ionic-covalent" bonding |Q. In fact, the CuO^ 
center with the valent {big^Aig — {bigeuYEu manifold represents a specific version of the "correlation" polaron, 
introduced by Goodenough and Zhou |32|] . 

The model allows to consistently explain many puzzling properties both of insulating and superconducting 
cuprates: the mid-infrared (MIR) region absorption bands the (pseudo)Jahn- Teller effect and related phe- 
nomena [Q, the spin properties |35|] . 

The presence of small polarons in semiconducting copper oxides has been detected with photoinduced infrared 
absorption measurements, infrared spectroscopy, and X-ray absorption fine structure techniques. 



One of the most exciting experimental evidences in favour of the model with the valent ^Aig — ^E^ multiplet 
is associated with observation in the doped cuprates of the MIR bands which polarisation features are compatible 
with those for ^Aig — ^E^ intra- multiplet dipolc transitions [ p9| . The corresponding transition energies observed for 
various cuprates are of the order of a few tenths of eV, that exhibits a typical energy scale for the valent multiplet. 

The Bu hole can be coupled with the big hole both antiferro- and ferromagnetically. This a rather simple 
consideration indicates clearly a neccessity to incorporate in valent multiplet both spin-singlet (&ige„)^£'„ and spin- 
triplet (6ige„)'^£'„, which energy could be even lower due to ferromagnetic big — e„ exchange. Indeed, the low- lying 
spin-triplet state for the two-hole CuO\~ center was revealed with the help of ^^'^^Cw NQR in La2CuQ^LiQ^04 



with singlet-triplet separation Ast = 0.13 eV |36 . Indirect manifestation of 02p7r, or e„ valent states were detected 
in Knight shift measurements by NMR for 123-YBaCuO system In connection with the valent ^Aig — ^E^ 
multiplet model for copper oxides one should note and comment the results of paper by Tjeng et al. | p^ , where 
the authors state that they "are able to unravel the different spin states in the single-particle excitation spectrum 
of antiferromagnetic CuO and show that the top of the valence band is of pure singlet character, which provides 
strong support for the existence and stability of Zhang- Rice singlets in high- Tc cuprates". However, in their 
photocmission work they made use of the (7^2^3/2(^3) resonance condition that allows to detect unambiguously 
only copper photohole states, hence they cannot see the purely oxygen photohole e„ states. 

It should be noted that the complicated ^Aig — -^'^Eu structure of the valent multiplet for the two-hole CuO|~ 
center has to be revealed in the photocmission spectra, all the more that the odd ^'^Eu terms play a principal role: 
namely these yield a nonzero contribution to ARPES for k = 0, or, in other words, at F point. In this connection 
one should note experimental measurements of the photocmission spectra in Sr2Cu02Cl2 and Ca2Cu02Cl2 
[ pO[ . All these clearly detect a nonzero photocurrent intensity in the BZ center, thus supporting the ^Aig — ^'^E^ 
structure of the ground state valent multiplet. 

1.2. Effective vibronic Hamiltonian and small (pseudo)Jahn- Teller polaron 

The orbital near-degeneracy leads to strong electron-lattice (vibronic) effects. For {^Aig,^ £'„) valent multiplet in 
the Cu04-center one has to consider the active Q-y nuclear displacement modes with 7 = aig, big, 629, e„. Effective 
vibronic Hamiltonian Hvibr is a sum of two contributions 

Hmbr — a EE + H AE ^ (1) 

that of vibronic coupling within the ^E^ term, and as well the vibronic coupling of the ^E^ and ^Aig terms, 
respectively. Vibronic Hamiltonian for isolated ^E^ term is familiar for a so called i? — 61 — 62 problem pa, and 
in terms of orbital operators (7 = aig, big, 629, e„) has a following form: 

ilEE = Vb,^ Vb,, Qbi, + Vb^, Qb^9 (2) 

with W6ig,62g being vibronic coupling parameters. Vibronic Hamiltonian Hae contains two terms 

HaE = Va,,K^Qa,,+Y.''e{QlV,^). (3) 

The first provides the connection between the A — E separation and symmetrical Qaig mode, while the second 
describes a linear vibronic coupling of the ^£'„ and ^Aig terms due to the active odd Qe^, Qel modes. One should 
note that for the CuO^ center one exists three types of the e„ modes. Vibronic Hamiltonian Hyib has to be added 
by elastic energy Uq for the CuO^ center 

^(Q)=E^' (4) 

i 

where summing runs over all normal displacements modes. A rather complete examination of the (pseudo)Jahn- 
Teller {}Aig,^ Eu)-ciig-big-b2g-eu problem was carried out in Ref. p^ . Depending on the relation between vibronic 
and elastic parameters one finds to be three situations: 

1. Weak vibronic coupling when Hmb results only in a rather small renormalization of bare elastic constants. 

2. The Jahn- Teller effect for the Eu term [E — bi — 62 problem) with emergence of two-well adiabatic potential 
for the rhombic modes with either big, or 629 symmetry. 

3. Strong pseudo- Jahn- Teller effect with vibronic mixing of (^^ig and ^E^) terms and emergence of four-well 
adiabatic potential. 

In practice one might observe various manifestations of the (pseudo) Jahn- Teller effect including the conven- 
tional static or dynamic effect, local structural instability, including the dipole one, spontaneous and induced local 
structural phase transitions with reconstruction of the multi-well adiabatic potential. 

The presence of the PJT centers, or small PJT polarons in semiconducting copper oxides has been detected, in 
particular, with photoinduced infrared absorption measurements, infrared spectroscopy, and X-ray absorption fine 



structure techniques. The photoinduced infrared absorption measurements on a number of insulating copper oxides 
have unambiguously revealed the formation of a localized electronic state accompanied by a localized structural 
distortion. These two aspects of the data have demonstrated the self-localized polaronic nature of the photo- injected 
carriers with a complex optical response spreading on a rather wide spectral range 0.1 1.0 eV. 

1.3. Large non-adiabatic hole polarons in the Cu02 layers of cuprates 

Non-zero hole transfer together with strong intercenter coupling due to the common oxygen favours the trans- 
formation of the small PJT polaron into a large non-adiabatic PJT polaron. Stabilization of a large polaron causes 
the hole to be shared over several CuO^ centers that could result in an optimal relaxation of the elastic energy and 
vibronic coupling all over the polaron volume. 

A model of large non-adiabatic hole polarons for the doped copper oxides was proposed by Bersuker and Goode- 
nough The large polaron could contain a hole cloud distributed over a significant number of the CuO^ centers. 
So, the authors |25| considered polarons containing 5 to 6 CuO^ centers. 

One of the model approaches capable to effectively describe such a polaron implies an extention of the conventional 
PJT problem with choice of a set of the symmetrized quasimolecular orbitals constructed from the valent hole states 
localized on the CuOa clusters involved in the PJT polaron (similar to the known LC AO-method), and a set of 
symmetrized displacements for the corresponding copper and oxygen atoms. 

The symmetrized quasimolecular orbitals can be build as linear combinations of the valent hole molecular orbitals 
tpjuir) for the CuOi cluster 

= E ^7^(1^ + r)V'7M(R + r)^o(R + r), (5) 

r 

where '!/'o(r) is ground state wave function of the CuO\^ cluster, FAf, 7/i are irreducible representations of the point 
symmetry group, C^ff are the symmetry coefficients, R is a radius-vector of the polaronic center of symmetry, 
and the sum runs over all the CuOa centers inside large polaron. For illustration we present the symmetrized 
quasimolecular orbitals for the large polaron containing 5 CuOa centers (see Fig. 2) with valent hig orbital: 

= x{o)M'j'^ - i{x(x) + x(-x) + x(y) + x(-y)}, 
= i{x(x) + x(-x) - x(y) - x(-y)}, 
= ^{x(x) - x("x) + x(y) - x(-y)}, 
= \{-x{^) + x(-x) + x(y) - x(-y)}, 

where |x) = IV'^ig^'big) = IV'oV'o) is the Zhang-Rice singlet wave function. 

Thus, we obtain a polaronic A\ — A\ — Bi — E manifold generated by the 6ig-state of the CuOi centers. This 
manifold incorporates two quasimolecular orbitals with s-symmetry, one with (ia;2-j^2-symmetry, and doublet with 
Px,2/-symmetry. The account for the hole transfer results in a bare splitting within the manifold. The orbital states 
with the same symmetry will mix forming the superpositions 

$rA/(R) -EE 4m^^7T(R + r)V^7M(R + r)^o(R + r). (6) 

One should note a specific bare quasimolecular orbitals for the A^-center "molecule" generated by local current 
states like e„-|- cx {e^x i icuy)- These manifest the A^-fold pattern of microscopic circulating currents resulting 
both in a state with nonzero magnetic moment and in a state with zero magnetic moment. In the latter case the 
time-reversal symmetry as well as rotational symmetry is broken but the product of the two is conserved. 

For many applications in scattering problems it is useful to introduce the form-factors 

/^^^Hk) = E4';';^^.Tw^*^ (7) 

r 

to describe the spatial distribution of the hole density inside a large polaron. 
Some form-factors for A^ = 4 large polaron (see Fig. 2) are 

/(^-)(k) = 2cos(^)cos(^), 



fi\fik) = -^V2sin(^)cos(^), fi\fik) = ^V2cos(^)sin(^), 
/^-^(k) = 2^sin(^)cos(^), /^^(k) = 2^^/2 cos(^) sin(^), 

/rx(k) = = 2 cos(^) cos(^), /-JJk) = /-J,(k) = 0. 

A set of the symmetrized quasimolecular orbitals $rM and symmetrized displacements Qtm for the copper and 
oxygen atoms inside large polaron form a basis for a rather standard, albeit complicated, vibronic PJT problem 
which solution gives an energy spectrum and appropriate wave functions '^aVn describing strongly correlated 
vibronic nature of the polaronic states. Here, aTn represent a set of quantum numbers which define a vibronic 
state. It should be noted, that only given the extremely simplifying assumptions this function could be written 
in a familiar Born-Oppenheimer form like ^ar-n = ^rMXaFniQ), where XaFniQ) is a vibrational function. In 
a wide sense, excitation spectrum of the large PJT polaron involves a whole vibronic spectrum originated from 
quasimolecular (7)r multiplet. It should be noted that the symmetry classification of polaronic hole states allows 
to elucidate some similarities with atomic system and partial waves; one might say about s— ,p— , ... like hole 
states. One should note an occurrence of the time-reversal symmetry breaking current states for the large PJT 
polaron. In addition, we must emphasize one more the specific role of the near-degeneracy for the valent manifold, 
and probable PJT effect, for the individual CuO^ center in formation of a large non-adiabatic polaron. 

Finally, large PJT-polaron in the lattice could be represented as a system of the CuO^ centers with a set of 
metastable states '^aVn specified by a binding energy and a life-time. Its nature implies strong charge fluctu- 
ations, so, it seems rather difficult to confine such a polaron within a single Cu02 layer. In other words, the 
large PJT polaronic nature of photo-hole implies its 31? structure with finite dimension, or correlation length, in 
the c-direction. Three-dimensional structure of the large PJT polaron could result in a rather strong fcj^, and 
consequently, photon-energy dependence of the ARPES intensity. Nevertheless, below we restrict ourselves, for 
simplicity, with the planar PJT polarons. 

In practice, for real systems the polaron will effectively couple with all the phonon modes which are active in 
the PJT effect Q. In other words, a large non-adiabatic PJT polaron may be considered as a bounded state of 
the large PJT-center and phonons. Coupling with phonons and spin system will result in an effective enlargement 
of polaron. The simplest way to account for this effect implies the introducing of a momentum-dependent cut-off 
factor to a polaronic form-factor. In addition, the phonon system will determine the relaxation dynamics of polaron 
states, and gives rise an effective dispersion to the hole spectral function. 

It should be noted that the quasiparticle behavior of polaron implies, as a rule, a single rather strongly bounded 
long-lived term. For this one can introduce an effective Hamiltonian which should include polaronic transport and 
coupling to phonon and spin lattice modes thus providing the coherent and incoherent part of the quasiparticle 
spectral function. It should be noted that the effective quasiparticle Hamiltonian could be look like as familiar Hub- 
bard, or t — J Hamiltonian. The remaining short-lived polaronic states will give rise to a rather wide structureless 
and dispersionless background in spectral region of valent manifold. 

Above we have considered the large lattice polaron and ignored the role of the antiferromagnetic background 
and appropriate spin fluctuations. Many authors have considered the formation of magnetic (spin) polarons p7| , 
moreover the assumption of spin-polaronic nature of the photo-hole, described in the framework of the extended 
t — J model is one of the most popular approaches to the interpretation of ARPES data in cuprates. So-called 
electron-hole asymmetric small polarons were introduced by J.E. Hirsch ||2^. In general, the polarons must be 
of complex spin-lattice hybrid type with a complicated spatial distribution of the electron and spin densities, and 
local structure distortions. 



2. Expression for photointensity 

Below we consider an expression for intensity of the photoemission with creation of immobile hole P JT polaron. 
An effective Hamiltonian for interaction with the electromagnetic field of frequency uj and polarisation e could be 
written within a polaronic manifold as follows 

H^nt = ^^AlrA/(k,e)cJ,>t^,, + i/.c., (8) 

TM k 

where k is a momentum of the final state of the photoelectron registered by the detector, c\^^ and ^p^v/^- ^-re creation 
operators for photoelectron and photohole, respectively. The matrix element is given by 

A^rM(k,e) = (^k(r)*[.l-'^|i/efl|^f )), (9) 

where 

HeR. = :^(p- A + A-p) 

Zmc 

is the interaction Hamiltonian with the electron momentum operator p and the vector potential A of the photon 
field; is the wave function for the ground state; ^Pp^ ^'^^ wave function for a TM state with one removed 
electron (one additional hole); 'i/'k(r) is the photoelectron wave function. It should be noted that expression (^) 
already implies a number of noticeable simplifications. 

Modelling the photoelectron wave function by a plane wave, we rewrite the expression for the matrix element 
(H) as follows 

A^rM(k,e) = ^ /(l^)(k)M^^(k, e), (10) 

where in the dipole approximation 

M^^(k,e) = (^,^(r)|(e.r)|e^'^'-). (11) 
Finally, the expression for the photoemission intensity may be transformed into 

/(k,w,e)cx ^ Alf^Mi(k,e)Alr2M2(k,e)AriAfi;r2Af2(w), (12) 
riAfi;r2Af2 

where the ground \g) and excited |e) states are the nonperturbed electron- vibrational states for the N-center cluster 
and vibronic ones for the hole PJT polaron, respectively. Emission spectral functions have a quite standard form 



^riAfi;r2Af2('^) = ^ XI ^ ^^'^^ {Ahl^M^a\9){9\hT2M2<y\e)S{uj + - Eg) = 

dte'^'ihhM.Amr.Aham- (13) 

Spectral functions contain a complete information about complex vibronic structure of the PJT polaron, and 
describe both the partial F-contributions at Fi ~ F2 and interference effects for different states with the same 
symmetry. These obey the sum rules 

^^riAfi;r2A/2(^) = "riA/i'5riAfi;r2A'/2- (14) 

Despite the extremely rough simplifying approximations Exp.(p^ displays very complex multi- level structure for 
photo-intensity with nontrivial polarization and k, lu dependence. 

Calculation of spectral functions AriA/i;r2A/2 (i^) for the large PJT polaron represents an extremely complex 
problem even at very strong simplifications p6[ |. For illustration, one might refer to a similar problem with spectral 
function which describes the line-shape of the optical A — E transition between the orbital singlet and orbital JT 
doublet 16). 

It should be noted, that, in a sense, the non-diagonal spectral functions Ar-^Mi;r2M2{^) describe the spectral 
weight transfer between Fi and F2 bands. In a whole, the PJT polaronic nature of the photo-hole provides a 
k-dependent ARPES spectral line-shape, in particular, with dispersive peak position. 



3. One-electron matrix element 



3.1. Copper contribution 

The copper atomic orbital with symmetry 7^ can be represented in the form 



where (7m) ^re coefficients specified by the symmetry requirements, Radi^) radial wave function, which we 
assume to be of simple Slater form 

2 AT r2 r r ^ 

R3d[r) = — W— 3 exp{~—}. 
81 V 15 a'l^/od iad 



Inserting these expressions to (|ll| ) and making use the familiar expansion for the plane wave |40 



e'kr^4^^ E ^''3L{kr)Yl,,{mLM{v) (15) 



L=OM=-L 



we rewrite (^11) as follows 



where we denote 



M(^")(k,e) = (d,^(r)|(e.r)|e*"-) = 



Dsik) = {R3dir)\r\Mkr)) ^ 62208\/| - 



5(l + 9a2fc2)5' 

4'^^)(e,k) = [Y^ X e']"'''* = J2 (-l)'^e_,CiTn,i^LM(k)ak(7/^). 

AI,q,m 

Here C]^-^_^ are the Clebsch-Gordan coefficients. The photoelectron energy dependence of quantities Di[k) and 
D^ik) given the Cu3d radial parameter = 0.35 A is shown in Fig. 3. 

3.2. Oxygen contribution 

The oxygen molecular orbital can be represented as a linear combination of atomic 02p functions centered at 
appropriate oxygen positions 

p^^ir) = J2 Cr(t)i?2p(|r - t|)ri„(r - t), (17) 

tm 

where C^J;^{t) are coefficients specified by the symmetry requirements, R2p(r) radial wave function. Inserting ([l7| ) 
to ( |lT|) and making substitution r' = r — t, we reduce (|ll] ) to 

M(k, e) = Y, Q.(t)(e • t)e**(i?2p(r')yi™(r')|e^'^"-') 

tm 

+ EC^;n(t)e'''*(i?2p(r')n„,(r')|(e.r')|e"^--'). (18) 

tm 

For convenience, one introduces two vectors with cyclic components 

Gm (k, e) = ^ C;„ (t) (e . t)^ , (k) = ^ ^™ (t)e^'^* . (19) 
t t 

Then Exp.(|l8|) could be rewritten in a more compact form 



M(k, e) = ^ G,„ (k, e) (i?2p (r)ri„, (r) | e"^"") + 

m 

J2 Zra{k){R2p{r)Y,ra{r)\{e ■ r)|e''''-). (20) 

m 

Making use the familiar expansion ( p^ ) for the plane wave we rewrite ( pO| ) as follows 

iV40)(k,e) = 2\/3^iS(fc)(G^^ • k)/fc+ 
{Ao(fc)(Z,, . e) - ^^(,) 3(e-k)(Z,,.k)^fc^(Z,,.e) ^ ^^^^ 

where we add the 7/i indexes label for the molecular orbital under consideration. Rather simple analytical expres- 
sions for the radial integrals Aq_2, B in ( pT| ) could be obtained if to make use the simplest Slater 02p radial wave 
function 

2v6 /qS ftp 



Then 



= (i?2p(r)|ji(fcr)) = 



^3 ■ 



Ao(fc) = (i?2p(r)|r|jo(fcr)) - 



3(l + 4a2fc2)^ 



(l + 4a2fc2)* 



512\/6a4 /a;fc2 

A2(fc) = (i?2p(r)Hj2(M>= ' 



(l + 4a2fc2)4 

The photoelectron energy dependence of quantities B{k), ylo(fc) and A2{k) given the 02p radial parameter ap = 0.52 
A is shown in Fig. 4. 

3.3. Expression for Aff,jg(k, e) 

Below we consider in detail the matrix element which specifies the contribution of the electron-removal process 
from the big orbital to form the Zhang-Rice singlet. 

For 7/x = hig at rather large photon energy Eph > 20 eV, but small binding energy £' < 1 2 eV) 



M,^;;"^(k,e) = 2i\l—{Di{k) + -D^{k)}{e,K, - CyKy), (23) 

I k ^ —— X . 

might obtain 



where Ha = ka/k, a = x,y,z. Making use the numerical values for the coefficients Cm^^\t) from Table 1 one 



^4^1? (^'^) = iy^{-3aB(fc)[K^e^cos(^) - Kj^e,, cos(^)] 



2[A^{k) + A2{k)][e, sin(^) - e, sin(^)] + 

&A2{k)[e^K.^ + eyK.y\[K.^s\n{^) ~ sin(^)]}. (24) 

It should be emphasized that the photocurrent intensity at the BZ center in the case of the even hig orbital as well 
as for any other even 7 orbital turns to zero. 

The bonding one-electron molecular hig orbital can be written as follows 

*f,i, (r) = db (r) sin Ob,^ -f pb (r) cos Ob,^ , 



where Oh^g is an angular covalcnt mixing parameter. Then the photocurrent intensity will be proportional to 

|(vl/,^^(k)|(e.r)|e*"-)p = [M^^^") p sin^ ^?,^^ + 
cos^ e,,g + sin^,,, cos9,,g{iA4y*Mi^J + A<"^(m(°V}. (25) 

3.4. Expression for Afe^(k, e) 

As it was mentioned above, the oxygen e„ states one might subdivide to a and tt orbitals. Due to a strong 
02p — 02p coupling they hybridize to form bonding and antibonding molecular orbitals p9[ |. The bonding 
one-electron molecular Cu orbital can be written as follows 

*e„^(r) =p(^^r) sin 6le„ + p'^^r) cos 6le„ , ^i = x,y, 

where 9e^ is a covalent mixing parameter. Then the photocurrent intensity will be proportional to 

J2\{^e.Ame ■ r)|e''^"-)P = J2 I^^^^^P sin^ + 

^ \Mjf^ P cos^ + sin 20,„ J2 H^'^H"^ ■ (26) 

Making use of the coefhcients from Table 1 one might obtain general expressions for matrix elements in the case 
of a states 7/i = CuX or 7/1 = e„j/: 

M('^)(k,e) = _V6^ai?(fc)e^K^sin(^)+ 

{{Ao{k) + A2{k))e^ - 3A2{k)K^{e ■ k)} cos(^) 

{^i = x,y). (27) 
Similarly, one might obtain general expressions for matrix elements in the case of tt states 7^ = e^x and 7/i — Cuj/: 

MW(k,e) = -V6^aB{k)eyK, sin(^)+ 
{(Ao(fc) + A2(fc))e, - 3A2(fc)«;,(e • cos(^), 

{{Ao{k) + A2{k))ey - 3A2{k)Ky{e ■ k)} cos(^). (28) 

4. Modelling the polaronic matrix element effects in ARPES spectroscopy 

Orbital big — e„ quasidegeneracy and polaronic nature of the photohole results in a complicated structure of 
the energy and momentum dependence of the photointensity. Below, we would like to present some examples of 
the straightforward model calculations of the matrix elements effects in ARPES. First of all we should address to 
the k-dependence of the single CuOi center contribution, the polaron form-factor effects, the photon polarization 
effects, and interference effects, caused by hig — e„ quasidegeneracy. 

4.1. Matrix element effects for isolated CuOi center 

The Figure 5 shows the contour-plots (the darker the color, the bigger the photointensity) for the quantities 
|Af^^(k)p which describe a partial one-center form-factor contribution to the photo-current intensity. The k- 
dependence of the polarization averaged |(^'fcjg(r)|(e • r)|e*'"')p for 61,^^ = — O.SStt is shown in Fig. 5a. Simply 
speaking, this is a contribution of the conventional model Zhang-Rice singlet. From left to right here we present 






the purely Cu3d, purely 02p, and the total contributions, respectively. One should note the complex k-dependence 
of the oxygen contribution as compared with the copper one. The k-dependence of the polarization and orbitally 
averaged |(5'e^^(r)|(e • r)|e''''')p contribution of the purely oxygen e„ electron-removal state for 0e„ = O.SOtt 
is shown in Fig. 5b. From left to right there are presented the n, a, and hybrid (interference) tt — a contributions, 
respectively. It should be noted that the tt and a partial contributions have a rather different k-dependence. 

To illustrate the photon polarization effects we present several examples of the angular (k) dependence of the 
photointensity for the "parallel" (e || k), and "perpendicular" (eJ- k) polarizations, respectively, calculated with 
the help of Exps. (23)-(24). The Figure 6a relates to the Cu3d partial contribution to the photointensity with 
creation of photo-hole in ZR-singiet state for the parallel (oc cos^ 20) , or perpendicular (oc sin^ 20) polarization, 
respectively. The Figure 6b relates to the corresponding 02p partial contribution. Here, one should note the more 
complex form of the angular dependence due to the essentially different structure of the respective matrix elements. 
However, the k modulus dependence of the polar plots in Fig. 6b is rather weak. So, in a whole, the polarization 
effects in both cases appear to be qualitatively, and even quantitatively similar. A comparative analysis of the 
Figures 5 and 6 indicates strong impact of the photon polarization effect on the final momentum dependence of 
the photocurrcnt intensity. 

4.2. Matrix element effects for isolated large polaron 

To illustrate an important role of the matrix element effects, we consider below a four-, and five-center model of 
the immobile large PJT polaron generated by the isolated valent ^Aig — ^'■^£'„ manifold of the CuO\~ center, that 
is assuming a localization of the photohole either in hig or e„ orbital on the CMO4 center. 

The Figure 7a, b shows the contour-plots for a number of quantities |/p'jj^''(k)p which represent a peculiar "fc- 
portrait" of the hole density within iV-center large polaron {N — 4, Fig. 7a, = 5, Fig. 7b), and describe a partial 
polaronic form-factor contribution to the photo-current intensity. Again we see complex and various momentum 
dependencies, reflecting both the hole symmetry and its distribution in large polaron. 

The Figure 8 shows the contour-plots for a number of quantities |A4rA/(k)p (F — Big,Eu) which describe 
a polarization and orbital averaged overall matrix element effect in partial F contribution to the photo-current 
intensity with one-center hole basis consisting of the big-, and e„-orbitals. The top figures present the total 
contribution, while below there are shown the partial big-, e„-, and interference big — e„-contributions, respectively. 
Interestingly, the interference term looks similarly in both cases. One should notice the nonzero contribution of 
the e„ states to the photointensity in the BZ center (F-point) for the E.^ type polaron. 

To illustrate the photon polarization effects we have calculated the angular (k) dependence of the partial pho- 
tointensity related to photo- hole creation in ^f^^ state for the "parallel", and "perpendicular" polarizations, 
respectively (Fig. 9). Interestingly, the dependence is qualitatively similar to the case of the ZR-singlet, at least for 
the high-symmetry directions. In other words, the polarization dependence alone could not distinguish the ZR- 
singlet among other terms with the relevant symmetry. The Figures 8 and 9 convincingly illustrate the role played 
by the matrix element effects, including the form-factor, orbital symmetry, photon polarization, and interference 
effects. As we see, for a quantitative comparison with experiment, it is necessary to add the ARPES amplitudes 
instead of intensities, taking into account the polarization and energy of the incident photons, and the direction of 
the photoemitted electrons. 

So, one might unambiguously say that the matrix element effects result in a complex k dependence of the 
photocurrent intensity which has to be taken into account when addressing such issues as Fermi surface. Indeed, 
the above model illustrations provide a wide choice of the "Fermi surface" -like behavior. Neglecting the matrix 
element effects results in erroneous conclusions concerning the electronic structure of the electron-removal states. 

5. Conclusions 

We had not for an object the detailed fitting of the experimental photoemission spectra, as we consider this 
problem in the meantime to be very complicated. At present, there is no generally accepted model for the large 
PJT polaron in copper oxides and appropriate ARPES spectral functions A{ll!), and it leads to uncertainties in 
quantitative interpretation of the experimental data. 

Nevertheless, we see that the low-energy ARPES spectra could be originated from the hole polaronic excitations 
which complex spectral shape is strongly affected by the soft lattice and spin fluctuations. These usually have to 
result in a complex ARPES spectral shape with a rather narrow purely electronic (coherent) peak and structureless 
(incoherent) background which describes the excitations accompanied by the emission and absorption of bosons 
(phonons, spin waves ). One should notice that the polaronic spectral response can spread over a wide energy 
range of about several tenths of eV. Experimental spectral shape of the intensity is, qualitatively, compatible with 
that expected for spectral response of the PJT polaron Perhaps, namely the polaron-like entity formation 

could explain the extremely narrow and intense peak lying below the Fermi energy and being the most intriguing 
feature of ARPES in all the high-Tc cuprates. We also could propose that many unusual features of the cuprate 
ARPES including the fc-dependence of spectral shape can be understood and described without any Fermi-surface, 
large or small. 



A qualitative comparison of experimental ARPES spectra and model calculations of the matrix clement effects 
allows to make a number of important conclusions: 

1. A model of the dispersionless valent ^Aig — ^'^E^ multiplet of the two-hole CuOf~ center and large (-/V=4, 
or 5) non-adiabatic PJT polaron is quite enough compatible with experimental ARPES data. Moreover, this is 
capable to describe some rather subtle spectral ARPES features, including the "remnant Fermi surface" effect 
without any reference to itinerant band-like states. It should be emphasized one more the principal role of ^'^E^ 
term determining the non-zero ARPES response in F point. 

2. The k-dependent spectral weight transfer for the large PJT polaron could be a natural origin of a "seeming" 
hole dispersion when to be described in a single-band model. At the same time, the polaronic transport described 
by the Hubbard-like model also could provide the natural explanation of the dispersion observed. Apparently, the 
straightforeward assignement of the photocurrent intensity maxima to that of the quasiparticle spectral density 
may result in erroneous conclusions. 

3. Interpretation of the photoemission spectra for strongly correlated oxides needs a caution and careful account 
for multi-band effects and matrix element effects, especially in what concerns the k-dependence of the spectral 
weight, and Fermi surface assignement. Moreover, the polaronic nature of the photo-hole gives rise to the problem 
of the "third dimension" of polaron, or to its correlation length in c-direction resulting in a number of important 
consequences concerning the photon energy dependence of photointensity. 

All these conclusions cast doubt on results of numerous papers with simplified interpretation of the low-energy 
ARPES data for Sr2Cu02Cl2 made in the framework of the various single-band versions of the t — J model 
[p|jic|], and aimed the interpretation of the "experimental quasiparticle dispersion law" supposedly determined 
from positions of the photocurrent intensity maxima. Unfortunately, the available experimental data do not allow 
to make so far the reliable conclusions on the nature of the low-energy ARPES feature. The ARPES data have 
to be considered with great care since it is most probable that they do not reflect straightforwardly the bulk 
density of states (DOS) of the quasiparticle excitations. In our opinion, a similar situation with the ARPES data 
interpretation occurs for many other strongly correlated oxides, including the superconducting cuprates, that forces 
to consider with caution many principal conclusions which are founded on the ARPES data. 

Concluding, one should be noted that the elaboration of adequate theory of the electronic spectra for the strongly 
correlated oxides needs at present not only solution of a number of complex theoretical problems, but the more 
perfect experimental data with complete spectral, angular, and polarization analysis. 

The research described in this publication was made possible in part by Award No.REC-005 of the U.S. Civilian 
Research & Development Foundation for the Independent States of the Former Soviet Union (CRDF). The authors 
acknowledge a partial support from the Russian Ministry of Education, grant # 97-0-7.3-130. 



Figure captions 

Fig.l. Schematic representation of the electron (hole) density distribution in the hybrid Cu3d02p big-, and 
purely oxygen eu((T, tt) molecular orbitals, assumed to be main components of the first electron-removal state in 
insulating copper oxides. 

Fig. 2. Schematic structure of the large polaron: (a) N — A, (b) N — 5. 

Fig. 3. Dependence of the CuSd-atomic radial parameters Di{E) and D^{E) on the photoelectron energy (a^ = 
0.35A). 

Fig. 4. Dependence of the 02p-atomic radial parameters B{E), Aq{E) and A2{E) on the photoelectron energy 
{ap = 0.52A). 

Fig. 5. Contour-plots for the /c-dependence of the partial one-center form-factor contribution to the photo-current 
intensity for the depolarized photons: a) jfj, — big, from left to right the Cuid, 02p, and the hybrid Cu3d02p 
{Ob^g = — O.Stt) contributions, respectively; b)7/i = from left to right the tt, a and hybrid cttt (0e„ = 0.47r) 
contributions, respectively. 

Fig. 6. Photon polarization effects. Angular {k) dependence of the ZR-singlet partial contribution to photointen- 
sity for the "parallel" (e || k), and "perpendicular" (e _L k) polarizations, respectively: a) CuZd partial contribution; 
b) 02p partial contribution. Numbers near curves indicate the k values. 

Fig. 7. Contour-plots for polaronic formfactors |/p]J^''(k)p which describe a partial polaronic form-factor contri- 
bution to the photo-current intensity, and represent a peculiar " /c-portrait" of the hole density within iV-center 
large polaron: a) iV = 4, b) = 5. 

Fig. 8. Contour-plots for quantities |A^rAf(k)P (F = Big,Eu) which describe a polarization and orbital averaged 
overall matrix element effect in partial F contribution to the photo-current intensity. The top figures present 
the total contribution, while below there are shown the partial big-, e„-, and interference big — e,i-contributions, 
respectively. 

Fig. 9. Photon polarization effects. Angular (k) dependence of the partial contribution to photointensity related 
to photo-hole creation in ^f^^ state for the " parallel" , and " perpendicular" polarizations, respectively. Different 
curves correspond the same k values as in Fig. 6b. 
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